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Abstract We study the Weyl-Wigner transform in the case of discrete vari¬ 
ables defined in a Hilbert space of finite prime-number dimensionality N. We 
define a family of Weyl-Wigner transforms as function of a phase parameter. 
We show that it is only for a specific value of the parameter that all the prop¬ 
erties we have examined have a parallel with the case of continuous variables 
defined in an infinite-dimensional Hilbert space. A geometrical interpretation 
is briefly discussed. 

Keywords Weyl-Wigner transform • Weyl-Wigner transform for discrete 
variables • Family of Weyl-Wigner transforms 


1 Introduction 

The Weyl-Wigner transform (WWT) was originally introduced to provide a 
phase-space representation of an operator defined in a continuous Hilbert space 
[TJ[5J[31I1] (see also [51[5] b When the operator in question is the density operator 
p, we speak of the Wigner function (WF) of the state. 

For a discrete Hilbert space of finite dimensionality, the WWT has been 
treated extensively in the literature. A representative list of references is given 
by [TlIRlini lTTnfTTlfriinTniTTlIT^fTHlITTini^ri 9 II 20 II 21 11^1^ : see also the references 
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contained therein. In Ref. m a definition was proposed which was appropri¬ 
ate for establishing a relation of WF with Kirkwood’s joint quasi-probability 
distribution and with the von Neumann model of measurement. 

Various authors (see, e.g., Refs. |10l[^l23) l have found a certain freedom in 
the election of the phase in the definition of the (finite dimensional) Wigner 
function. In the present paper we show that the definition for the discrete 
case given in Ref. m can also be extended to define a family of WWTs as 
function of a phase parameter c. We study this freedom in detail and show 
that one particular case, c = 0, reproduces the definition given in Ref. |22) . 
while another member of the family, c = —1/2, brings the various expressions 
for the discrete case to a form similar to that for the continuous case. 

The paper is organized as follows. In Sec. [5] we define the family of WWTs 
for the discrete case, describe some of its properties, and analyze the special 
role played by the choice c = —1/2 in the structure of the discrete WWT when 
it is compared with the continuous case. The approach we adopt in this pre¬ 
sentation provides a somewhat different view of the intimate role (discovered 
by Grossmann and Royer [MHH]) of parity in the phase-space formulation of 
QM. In Sec. [3] we present an alternative form of the WWT for the continuous 
case, which exhibits its relation to the discrete one from another angle. In Sec. 
|4] we briefly discuss a geometrical interpretation of our approach and com¬ 
pare it with that studied by other authors. Sec. [S] contains some concluding 
remarks. For coherence and completeness, and in order not to interrupt the 
main flow of the presentation, we relegate to the appendices some proofs and 
discussions of properties of various concepts used in the text. 


2 Discrete Weyl-Wigner Transform 

2.1 Definition and properties 

The possibility of defining WWT for a Hilbert space of finite dimensionality 
has been studied by many authors [TieiaiiniiiiiiiiiiiiiiiiniiisiiiiiiiHiiiijiini 
Here we extend the definition of Ref. [32] and define a family of 
WWTs: using the Schwinger operators Z and X (see Ref. [33], where these 
operators are designated as U and V) for a Hilbert space of a prime-number 
dimensionality N (summarized in App. we define, for an operator A, and 
as function of the parameter c, the family 

.. ( N-lN-l 

f 6=0 fe=l 

N-l . 'j 

+ Tt ^ } ■ (1) 

We have defined a discrete phase space in which the “coordinate-like” and 
“momentum-like” variables are denoted hy q^p = 0,1, ■■■, N—1. The operators 
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[N{N — 1) in number; 6 = 0,1; k = 1, ■ ■ ■, N — 1], together 

with the N operators {k = 0, ■ ■ ■, N — 1) which appear in Eq. (ITJ, form a 
complete set of N'^ operators. Notice that this set of operators is the finite¬ 
dimensional analog of the continuous complete set eyip[i{uq -I- vp)]. If N" is a 
prime number, the integers 0,1,... ,iV — 1 form an algebraic field analogous 
to that of the real numbers in the continuous case. The N = 2 case requires 
a special treatment. We note that the exponents of a; = exp(27ri/iV), which 
will appear frequently in our analysis, always belong to the Mod[7V] algebra. 
We remark that, when iV is a prime number, there is a simple expression for 
the mutually unbiased bases (MUB) [27) . as summarized in App.lBlfsee. 

e.g.. Refs. [IZlIlIlliH])- 

The definition o can be given an alternative expression, whose interest 
lies in the fact that it makes its manipulation easier and, also, suggests an 
interesting geometrical interpretation: the expression (HD can be written in 
terms of MUB as 


N-l N-l N-l 

Wf{q,p) = ^ E 

b—0 m—0 


N-IN-1 

E E (rn;b A m;b'j - Tr(i) , (2) 


where the reference basis (eigenstates of Z) has been denoted, for convenience, 
as (i, and we have defined 

( q, for 6 = 0 . 

In Eq. ([2D we can do the sum over k, so Eq. (|T|) can be given the alternative 
form 


Wfiq.p) = E ^ mW( 6);6) -Tr(i) , 


N-l 


6=0 


= Tr(i4(;)) , 

where we have defined the Hermitean operator 


(4) 

(5) 


p(c) 

Q^P 


N-l 

E |MW(6);6)(MM(6);6|-i, 


b=6 


( 6 ) 


(I being the unit operator), which we now analyze. For a given pair of phase- 
space variables q,p, the expression Mq'^plb) may be considered as defining 
“points” in a b-m plane: 6 is along the horizontal axis and takes the values 
6 = 0,0,1,---,A^ — 1, which denote the -|- 1 bases; m is along the vertical 
axis and takes the values m = 0, — 1, which denote the N states for 

each basis. Thus, for fixed q^p, this set of points may be taken as defining a 
“line” in the b-m plane (see Fig. |T]for an example): we thus refer to the set 
Mq‘^l{b) (6 = 0, 0, • • •, — 1) as a line, and to the operator Pq'^p of Eq. ([BD as 

a line operator. 
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Our first aim is to obtain the matrix elements of the line operator. We split 

p{c) ^ rp{c)y , p/r 

q,v \ q.p / ' q.p 1 \'/ 


Pq^^l as 


where 


N-1 

= E |mW(6); 6)(mW(6); b\-I, (8) 

b=0 

P; = |g)(g|, 6 = 0. (9) 

We refer to as the “amputated” line operator and consider it first. 

A diagonal matrix element with respect to the reference basis of any one 
of the summands in the first term of Eq. ([81) equals 1/iV, according to the 
definition of MUB, as can also be verified from the explicit form of the MUB 
states given in Eq. (HSl) . Thus the amputated line operator {Pjpl)' has zero 
diagonal elements, i.e., 

H(Pip^))'|n) = 0, (10) 

because of the subtraction of I. 

We now consider a non-diagonal matrix element for an arbitrary 6. From 
Eq. (H51) we hnd 

(n|M,(;)](6);6)(6;M,(;)](6)|ri') = 

If we compare this matrix element at 6 with one at b' p b, we see that the two 
will be equal iff 

-= (g + c)Mod[iV], (12) 

a condition to be referred to as the coherence requirement. When we add up 
terms like (HU for 6 = 0, • • •, — 1 to construct the n p n' matrix element of 

(|8|). only matrix elements with n + n' = 2g + 2c + 1 will add up coherently; 
other matrix elements will cancel each other as they run over the roots of unity. 
Substituting, in Eq. (ITTl) . n + n' = 2q-\-2c-\-l and Mg‘lp{b) = —p + b{q + c) 
from Eq. we obtain, for the matrix elements {n\(Pqp'^)'\n'), N equal terms, 
with the result 


(n|(p(^))'|n') = W.2,+2c+i n p n'. 


For the matrix elements of P^p we have 

{n\Pqp\n') = SnqSn'q ■ 

The matrix elements of the full line operator are then 


{n\P^;'>\n') = |^-+-'.29+2c+i w 


{n—n')p 


n p n' 
n = n' 


(13) 


(14) 


( 15 ) 
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which can also be written, Vrijn', as 

(n|pM|n') = - Snn' ,+e+i + W,2,+2C+1 (16) 

The structure of the matrix {n\Pqp\n') of Eqs. (ITSl) . (fTBl) . is illustrated for 
iV = 5 in Tables [T] and [21 for the values c = —1/2 and c = 0 for the phase 
parameter, respectively. For c = —1/2, the non-zero matrix elements for fixed 
q appear precisely on the various secondary diagonals; notice that, for a given 
g, the diagonal matrix element SnqSn'q lies on the same secondary diagonal as 
the off-diagonal matrix elements; this is not the case for c ^ —1/2. This shows 
the particular role played by c = —1/2. We thus see that the phase freedom 
represented by the parameter c has important consequences in the structure of 
the matrix of the line operator. We remark that this phase freedom can only 
be there for the discrete case; in the continuous case this is not allowed, if we 
insist that the resulting function be continuous. It is thus the value c = —1/2 
for the phase parameter in the discrete case that shows a similarity with the 
continuous case. 


^ n / n' ^ 

0 

1 

2 

3 

4 

0 

^nO^Ti'O 


la = l| 


II 

to 

1 




9 = 2 


2 

|g = i| 


^n2^7i'2 


II 

CO 

3 


II 

to 


<5n3<Sn'3 


4 

9 = 2 


CO 

II 


^7T4^n'4 


Table 1 Structure of the matrix {n\Pq'^\n') of Eqs. 11511 . 11611 . for N = 5 and c = —1/2. 
The non-zero matrix elements for fixed q (which obey n -|- n' = 2q) appear on the various 
secondary diagonals. For a given q, the diagonal matrix element SnqS^/q (i.e., n = n' = q) 
appears precisely on the same secondary diagonal as the off-diagonal matrix elements; this 
is emphasized for the particular case 9 = 1, where the corresponding matrix elements have 
been highlighted and a box has been drawn around each one of them. 


An important property of the operators Pjpp is that they form a complete 
orthogonal set of operators. For any value of c: 
i) They fulfill the orthogonality relation 


—p('=) 

jq [ <i,P q',p' 


— ^q,q' ^p,p' 


(17) 


This statement is proved in Add. (Cl employing geometrical arguments, 
ii) They satisfy the closure relation 


1 ^-1 

( 18 ) 

q,p=0 

This follows from Fq. m- for n ^ n', the sum over p gives zero; for n = n' 
we have J^qp ^nq = N. 
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4- n / n' ^ 

0 

1 

2 

3 

4 

0 

^nO^n'O 

O 

II 


|g = l| 

q = 4 

1 

q = 0 




q = 4 

q = 2 

2 


|i'? = i| 


&n2&n'2 

q = 2 


3 

|9 = 1| 

II 

q = 2 

SnsSn'S 

II 

CO 

4 

II 

q = 2 


>0 

II 

CO 



Table 2 Structure of the matrix {n\Pqp^\n') of Eqs. lllSIl . I I16I1 . for N = 5 and c = 0. The 
non-zero, off-diagonal matrix elements for fixed q (which obey n -|- n' = 2q 1) appear on 
the various secondary diagonals. For a given q, the diagonal matrix element SugS^'q (i-e., 
n = n' = q) does not lie on the same secondary diagonal as the off-diagonal matrix elements; 
this is emphasized for the particular case q = 1, where the various matrix elements have 
been highlighted and and a box has been drawn around each one of them. 


iii) An N x N matrix A can thus be written as a linear combination of the 


n(c)5 


i.e., 


^ N-l 

^ = V S Tr(iAyS 

(19) 

q,p=o 


N-l 


q,p=o 

(20) 


where we have used Eq. 

From the above discussion we find that the WWT of Eqs. o and ® 
possesses the following properties. 

1) For A = p, the WE is normalized as 

N-l 

- W^^\q,p) = l. (21) 

p,q=0 

2) The WWT of a Hermitean operator A is real, i.e., 

W^f{q,p) = [W^f{q,p)r, for At = A (22) 

This follows immediately from the Hermiticity of the operators Pq'^p. 

3) The WWT of the unit operator is 1, i.e., 

Wf(g,p) = l. ( 23 ) 

4) The WWTs of the operators A and B fulfill the so-called “product 
formula” (see also Ref. [TO], Eq. (15)) 

AT-l 

E w[^\q,p)wl^^\q,p)=TT{AB) . 

q,p=0 


(24) 
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5) The WF W~‘^\q,p) satisfies the marginality property, written in terms 
of the projector = \mb){mb\, 

^ N-l 

Tr(pP^b) = {m,b\p\m,b) = ^ ^ ’ (25) 

q,p=0 

where we recall that Mq'lp{b) is defined in Eq. Q. 

Eq. ([251) states that the probability to find the system in the state m of 
the basis b (of our set of TV + 1 MUBs) is 1 /N times the sum of the WF over 
Mq'^pib) = m. Its RHS can be considered as defining the Radon transform of 
the WF wj:‘^\q,p) (see, e.g.. Refs. [51151 1^150] 1. Two particular cases of the 
above marginality property are 

(9o|pko) = , (26) 

P 

{Po\p\Po) = ^t,‘'\9,Po) , (27) 

9 

which are the standard marginality relations. 


2.2 The particular case c = —1/2 

We aready remarked that the choice c = —1/2 for the phase introduced in the 
original definition m constitutes an important particular case of the general 
formalism. We now examine this choice in relation with another property of 
the WF. 

In this case, the matrix elements (ITSl) . (IT51) of the operator Pq'^p reduce to 
the simpler form (we remove, for brevity, the c = —1/2 index in the various 
expressions) 

(n|P,p|n') = Vn,n'. (28) 

For this choice, the operator Pq p can be related to the parity operator Ilq^p 
around the point {q,p) as [2411251131] 

Pq,p = Hq.p ■ (29) 

This is quite evident, e.g., when q = p = 0, from Eq. (1551) . Indeed, 

{n\Po.o\n') = Sn+n ',0 (30) 

N-l 

and hence Pop = \n){-n\ , (31) 

n—0 

which is the parity operator around the point {q,p) = (0,0). More generally, 
the operator 

nq,p = x^zppopZ-px-\ 


(32) 
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which can be seen to have the same n^n' matrix elements as Pq,p: can be 
verified to correspond to the parity operator IJq^p around the point (q, p). One 


can also easily verify that {nq,p)^ = I, further verifying its shifted parity 


character. 

This property, which we showed for c = —1/2, has a parallel in the WWT 
for the case of a Hilbert space for continuous variables UMrnn- 

The discussion given in the next section gives a better understanding why 
the choice c = —1/2 in the family of discrete WTTs allows us to achieve, 
in a natural way, a structure for the WWT which is similar to that for the 
continuous case. 

3 Relation with the WWT for continuous variables 

We first re-write the definition m of the WWT for discrete variables in terms 
of position-like and momentum-like operators using (HU) and (PI) , as 



gi^fc[-p+fc(g+c)] 





(34) 


We shall use this expression below. 

We now turn to the case of continuous variables. The WWT can be ex¬ 
pressed as the inverse Fourier transform of the characteristic function of the 
operator as [dllBlI^ 



(35) 


where we have used units in which q and p are dimensionless, and h = 1. 

It is shown in App. |D]that the definition can be written in the alter¬ 
native form 



s[—p-l-cot 9 (<J— 5 )] 


ir sin 


(36) 


xe' 


Notice that, having expressed the WWT in terms of the operator e“®P e^cote q 
as in Eq. (pil . [instead of the operator used in Eq. (IHKl) ]. there 

appears naturally a shift —1/2 in the coordinate q. 
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We now compare Eq. (IMl) for discrete variables with Eq. (IMll for the con¬ 
tinuous case. Notice that the first term in Eq. (IMl) has the structure of Eq. (IHf)!) 
for the continuous case when 0 ^ 0, ii we make the correspondences shown in 
Table [31 For 0 = 0, we notice that 

lim I (e“*^ gi cot 9 9 ^ ® gir sin 9 [-p-|-cot 9 (g- i )] | _ ^-irq^irq ^ 1 - 0 ^^ 

which has the structure of the second term of Eq. (1341) , identifying k with r. 
However, the fact that the —1/2 shift in q goes away for 0 = 0 is of no con- 


Eq. UMI 

Eq. [13611 


Q — ip 0'i cot 9 q 

b 

cot 0 

k 

r sin 0 

c 

-1/2 


Table 3 Correspondence of the discrete and continuous WWT’s, for 6^0. 


sequence in the continuous-variable case, as it occurs in a set of zero measure 
and thus does not affect the integral (36). This is consistent with our previous 
discussion in relation with Tables 1 and 2, where we noted that the phase 
freedom that we have in the discrete case is not there in the continuous case. 

The above discussion gives a better understanding why the choice c = 
— 1/2 for the discrete case gives a WWT having similar properties as for the 
continuous case: in point of fact, in the latter, continuous case, a shift —1/2 
in q occurs in a natural way when we express the operator in terms 

of the product of operators e“®^ gjcot9 q g^g g-.^g ^ similar structure as 

in the former, discrete case. 


4 Geometrical interpretation of the “phase freedom” and its 
relation with other approaches 

As we mentioned in Sec. O the considerations briefly mentioned there lead us 
naturally to a geometrical interpretation of our analysis (see, e.g., Ref. [3T]1. 

Notice that our approach is based on a “dual affine-plane geometry” (DAPG) 
(where the points represent states), in contrast to the “afRne-plane geometry” 
(APG) (where the lines represent states) used, e.g., in Refs. (TU], and [20] ■ In 
Ref. [20| both approaches are discussed. 

i) In the present (DAPG) approach one has N{N + 1) points and lines. 
The points represent MUB projectors, and the lines, Wigner operators. The 
lines defined by m = Mq'^l{b) depend on the freedom to choose the phase 
parameter c. This dependence is illustrated in Fig. |T] for the particular cases 
c = 0 and c = —1/2, the latter having played a particular role in our earlier 
discussion. The matrix elements of the line operator, Eqs. (HU), (HOj), reduce. 
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Fig. 1 Illustration of the line m — in the b — m plane, for the dimensionality N — b and 

the particular pair of “phase-space” values q — 2 and p — 1. The dots correspond to the choice 
c = 0 of the phase parameter and the crosses to c = —1/2. 


when the choice c = —1/2 is adopted, to the “standard Wigner operator” (1^ . 
a relatively simple expression which, for convenience, we reproduce here 

(n|P,p|n') = Vn,n'. (38) 

In contrast, the choice c ^ —1/2 leads to the general ’’non-standard”, more 
complicated, Wigner operator of Eqs. ca), lUl). Notice, however, that the line 
operators are orthogonal Vc, as shown in Eq. (na. 

We thus see that the phase ambiguity arises via the ambiguity in defin¬ 
ing the line equation Mqp{b), giving “standard” or “non-standard” Wigner 
operators. 

ii) In the (APG) approach of Ref. [TU] one contemplates N'^ points and 
N{N + 1) lines. One has a freedom in defining the point operators Aa, a = 
{q,p), SO that they are “standard” or ’’non-standard” Wigner operators. The 
latter are similar to the “non-standard” Winger operators we get in the DAPG 
approach when cA “1/2. The lines give MUB projectors. 

We thus see that the phase ambiguity arises via the ambiguity in defining 
the point operators, such that they give “standard” or “non-standard” Wigner 
operators. 

In addition, a particular freedom mentioned in Ref. m, p. 11, is a linear 
transformation of the phase space coordinates, leaving their basic Eqs. (Al’)- 
(A3’) intact. 

Ref. [20] (Sec. 3.6) finds, too, a freedom in the election of a phase in the 
definition of Wigner’s function for the discrete case. The authors use an APG 
approach, so we can apply arguments similar to those used above for the APG 
approach of Ref. m- In addition. Ref. m describes a phase freedom which. 
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to the best of our knowledge, is not the same as the one found in the present 
paper. The authors discuss the freedom in the choice of the “rotation operator” 
in the plane: this rotation transforms among the bases. In our approach, the 
various bases are “vertical lines” (see Fig.[T]) in the b-m plane; we have chosen 
for them the form given in Eq. (I45L which is kept fixed in all our analysis, 
without further changes in its structure. Also, Ref. [50] finds a relation with 
parity which seems to be independent of the choice of their phase (j>. In contrast, 
in our case, similar relations hold only for the choice c = —1/2. 

5 Conclusions 

We considered the mapping of Quantum Mechanics from Hilbert space to 
phase space, which is arguably the most intuitive form of the three Quantum 
Mechanics formulations: Hilbert space, summation over paths and phase space. 
We demonstrated that for a given finite-dimensional Hilbert space there are 
several phase-space forms possible which differ in what we termed a phase 
parameter. We showed that one particular value of the phase parameter brings 
the various expressions for the discrete case to a form similar to that for the 
continuous case. In particular, our approach provides a novel view of the role 
of parity in such mappings: the “most natural” choice for the phase parameter 
gives rise to parity, in a way that parallels the situation for the case of a 
Hilbert space for continuous variables. We briefly compared the phase freedom 
discussed in this paper with that found by other authors. 

Acknowledgements One of the authors (PAM) acknowledges supposrt by DGAPA, under 
Contract No. IN109014. He also acknowledges the kind hospitality of the Physics Department 
of the Technion, where this investigation was initiated. 

A Schwinger operators 

We consider an Af-dimensional Hilbert space spanned by N orthonormal states |g), with 
q = 0,1, • • •, (Af — 1), and subject to the periodic condition |q-|-Af} = |q}; they are designated 
as the “reference basis”, or “computational basis” of the space. We introduce the unitary 
operators X and zm, which fulfill the periodic condition = I (I being the unit 

operator) and are defined by their action on the states of the reference basis as 


Z|q} = a;«|q), u = , 

X\q) = |q-|- 1), 


(39) 

(40) 


leading to the commutation relation 


(41) 


ZX = uiXZ. 


The two operators Z and X form a complete algebraic set |2fi| , so that any operator defined 
in our A^-dimensional Hilbert space can be written as a function of Z and X. 

We introduce the Hermitean operators p and q through the equations 



X = UJ-P = , 

Z = uj^ = . 


(42) 

(43) 
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As X performs translations in the variable q and Z in the variable p, we designate p and q 
as “momentum-like” and “position-like” operators, respectively. 

What we defined as the reference basis can thus be considered as the “position basis”. 
With (gll) and definitions ilia . Il43l l. the commutator of q and p reduces, in the continuous 
limit 128)1320331 . to the standard one, [q-,p] = i. 


B MUB 

The operators XZ^, b = 0^ • ■ ■ N — 1 define N of the TV + 1 MUB, [see Eqs. 1144 l l. (I45I I below], 
while the operator Z defines the reference basis. The operator XZ^ possesses TV eigenvectors, 
denoted by \m^b) (see Ref. |17| and Eqs. (10), (11) of Ref. |21|) 

XZ^\m,b) = Lj'^\m;b); fe, m = 0,1, • • •, TV - 1 , (44) 

N-l 

I™;'’) = (45) 

^ q = 0 

Here, Iq-) (q = 0, • • •, TV — 1) denote the TV states of the reference basis. The states with 6 = 0 
are eigenstates of p 


|m; 0) 


_L y 

\p = —m = (TV — m)Mod[TV]). 


(46) 

(47) 


C Orthogonality of the operators Eq. (I17|) 


We show the orthogonality, as given by Eq. 03, of the operators Pq^p defined in Eq. IS}. 
We first note the following two features of the lines described in the text. 

1. We assume the phase c to be fixed; then a line is defined by the pair of numbers q,p. 
Lines Pq^p, Pqf^pf may be either identical, which means that p = p' and q = q\ or distinct, 
meaning that either 

i) 9 ^ g', P = p', or 

ii) q = q\p^ p', or 

iii) q ^ q' and p 7 ^ p'. 

If two lines are distinct, then they have one point in common, i.e. there exists (one) 6 
wherein the two lines intersect. Every line has one common ’’point” with every other line. 
If two lines have two points in common, they are identical. 

Let line 1 be given by {q,p) and line 2 by (q',p'). We have the possibilities: 
i, iii) The lines differ if q 7 ^ q^ Here both cases p = p' or p 7 ^ p' imply 6 = (6 = 0 for 

the first case) - a unique value in either case. 

ii) If the lines differ via p 7 ^ p' but q = q', the common point is at 6 = 0. 

The uniqueness of the solutions imply that distinct lines cannot have more than one 
point in common. 

2. All the “points” are MUB projectors , i.e. 

p(")(g,p; 6 ) = |MM(6);6)(M.';4(6);6| . (48) 

For 6 7 ^ 6 ' (including 6 = 0) we have 

[pW(q,p;6)p('=)(g,p;6')] = ^ ■ 


Tr 


(49) 
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The actual proof of orthogonality involves computing the LHS of Eq. (I17l> . giving 


N-l 

Tr = E Ti- + E Tr &')] 

b=6 

N-l N-l 

-ETr[p( ^\q,p-,b)^ - Tr |^p(‘')(i}',p';6)j + Tr I, (50) 

6=0 6=6 

= A + B + C + D + E. (51) 

Calculation of A . The first term, A, involves a sum over N 1 traces of the product of 
two projectors, both in the same basis b. 

i) li q = q',p = p', the lines are identical: then the sum reduces to V + 1 traces of 
projectors, each giving 1, the result thus being N 1. 

ii) If the lines are distinct, they have only one common point at b; we thus have the trace 
of a projector, giving 1. All other N terms involve the trace of a product of two orthogonal 
projectors and do not contribute. Thus 

A = ^TrP(")(g,p;6)p(")(g',p';6) = {N + l)Sg^,5pp, + 1 ■ {1 - = 1 + AfVV • 

b 

(52) 

Calculation of B. The second term involves traces of projectors of distinct bases. Since these 
are MUB projectors, each term gives 1/N. There are N(N + 1) terms in the sum, giving 
N{N = N ^1: 


B= = N+ 1. (53) 

bj^b' 

Calculation of C and D. The third and fourth terms, C and D, involve traces of projectors 
each multiplied by unity. There are A + 1 terms, each giving 1. We thus find 

C = -y^TrP(‘=)(q,p;fe) = D = ,p';b') = -{N + 1). (54) 

b b 

Calculation of E. The last term, E, involves the trace of unity, giving N. 

Adding up the five terms, we finally obtain: 

Tr ■ (55) 


D Derivation of Eq. (|36p 

The definition of WWT, Eq. II35II . is based on the notion that the Weyl operators [I] 

U = (56) 

form a complete and orthogonal operator basis t5ll6l[30l. For the purpose of comparing the 
continuous and discrete cases in Sec. O we introduce the alternative set of complete and 
orthogonal operators 

V = e»“'3 , (57) 


and express the WWT in terms of them. 

Making the change of variables u' = —v and relabelling v' again as u, we obtain, from 
Eq. (I35II 


1 POO POO 

W,iq,p^ = -TrLJ_ 


— oo </ —oo 




(58) 
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Mann, Mello, Revzen 


We use polar coordinates 


and write 


u = r cos 0 , V = r sin 0 


roo r2‘K 

W^{q,p) =Tr rdr ^e-iraine{-23+cot S q) 

7o 7o 

^ gir sin 6(—p+cot B q) 

= Tt rdr fj d9 A « «] 

sin 6(—p+cot 6 q) 

Recalling the BCH identity 


and choosing ^ = — 1, a = cot0, we can write 

gi( —p+cot 8 q) _ g^ cot 0g —ip gi cot B q 


,i( —p + cot B q) 


— ( c. — '^P (J- cot B q 


-\ —r sin 6 


so II61II takes the form of Eq. II36I I given in the text. 


(59) 

(60) 

(61) 

(62) 

(63) 

(64) 
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